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ABSTRACT 


We report on NuSTAR and Swift observations of a soft state of the neutron star low- 
mass X-ray binary GS 1826-24, eommonly known as the “eloeked” burster. The tran¬ 
sition to the soft state was reeorded in 2014 June through an inerease of the 2-20 keV 
source intensity measured by MAXI, simultaneous with a decrease of the 15-50 keV 
intensity measured by Swift[QPIA. The episode lasted approximately two months, after 
which the source returned to its usual hard state. We analyze the broad-band spectrum 
measured by Swi/t/XRT and NuSTAR, and estimate the accretion rate during the soft 
episode to be « 13%mEdd> within the range of previous observations. However, the 
best fit spectral model, adopting the double-Comptonization used previously, exhibits 
significantly softer components. We detect seven type-I X-ray bursts, all significantly 
weaker (and with shorter rise and decay times) than observed previously. The burst 
profiles and recurrence times vary significantly, ruling out the regular bursts that are 
typical for this source. One burst exhibited photospheric radius expansion, and we 
estimate the source distance as (5.7 + 0.2) kpc, where fb parameterizes the possi¬ 
ble anisotropy of the burst emission. Interpreting the soft state as a transition from an 
optically thin inner flow to an optically thick flow passing through a boundary layer, 
as is commonly observed in similar systems, is contradicted by the lower optical depth 
measured for the double-Comptonization model we find for this soft state. The effect 
of a change in disk geometry on the burst behavior remains unclear. 

Subject headings: accretion - X-rays: binaries: close - stars: neutron - X-rays: bursts - 
X-rays: individual: GS 1826-24 
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1. Introduction 


Type I X-ray bursts arise from unstable thermonuelear burnin g on the surfaee of aeereting 


neutr on stars (NSs) in low-mass X-ray binaries (LMXBs; see, e.g., 


Lewin. van Paradijs & Taam 


19931 for a review). While variations in burst properties from souree to souree, and with time, 


are explained by ehanges in the aeeretion rate and the fuel eomposition at ignition, a detailed 
physieal understanding of most X-ray burster systems is still laeking. For example, there is no 
explanation for the deerease in burst rates (leading apparently to a transition to stable burning) 
that oeeurs at aeeretion rates about a faetor o f ten below the theoretieally expeeted value (see, e.g.. 


Cornelisse et al.l 


2003 


Galloway et al 


20081). The details of the relationship between aeeretion 


rate, burning physies, burst morphology, and burst reeurrenee times are eomplex and still not 
understood. 

For the majority of burst sourees, that aeerete a mix of hydrogen and helium from 
their eompanion, a general pieture of bursting behavior arises with four burning regimes 


marked by inereasing loeal a eeretion rate (m) per NS unit area (see, e.g.. 


Fuiimoto et al. 


Strohmaver & Bildsten 


20061 for details). As identified in 


Fujimoto et al. 


1981 


198lh . ease 3 burning 


oeeurs at low (< 0.01 mEdij) aeeretion rates, and arises from unstable hydrogen ignition in a mixed 
H/He environment. No example of this type of burning has eonfidently been observed. At higher 
aeeretion rates, eorresponding to ease 2, steady burning of hydrogen eommenees, while helium 
burning is still unstable. However, m is low enough that the aeereted hydrogen is exhausted at the 
base of the fuel layer by the time unstable helium ignition is triggered, so ease 2 bursts should 
oeeur in a He-rieh environment, and with relatively long reeurrenee times. The resulting burst 
light eurves exhibit short (< 1 s) rises and tails (< 10 s), with high peak luminosities typieally 
exeeeding the Eddington limit, whieh generates photospherie radius-expansion. Henee, the ratio. 


^Defined as the mass aeeretion rate, 8.8 x 10"* g em ^ s \ eorresponding to the Eddington lumi¬ 


nosity in a 1 AMq NS frame. 
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defined as the a value, of the persistent fluenee between bursts to the burst fluenee often exeeeds 
100. At higher aeeretion rates (above a few pereent of mEdd) the burst reeurrenee time beeomes 
short enough that hydrogen remains in the base of the fuel layer at ignition, and these ease 1 bursts 
exhibit long profiles eharacteristic of jS-deeay-mediated hot-CNO burning, and rp-pwcess tails. 
The a values are eonsistently lower than in ease 2. Finally, at aeeretion rates > mndd, He-buming 
should stabilise, and no further bursts are expeeted. 


GS 1826-24 (aka Ginga 1826-238, the “eloeked” or “textbook” burster; see 


Ubertini et al. 


1999h demonstrates the elosest agreement with theoretieal model predietions among the over 
100-known thermonuelear burst souree^. It has exhibited regular bursting behavior with highly 
eonsistent properties from burst to burst over the 30 years since its dis covery as a new transi ent 


(Tanaka et al. 1989). Indeed, using RXTE observations of 24 bursts. 


Galloway et al. 


(2004) 


measured a relationship between persistent X-ray flux and burst recurrence time: the latter 
decreases almost linearly as the accretion rate increases. This implies that the accreted mass 
between two bursts is each time completely burned during a burst, and is approximately the same 
even as the accretion rate changes. The burst light curves and properties of GS 1826 -24 have also 


been shown to be in good agr eement with the pr edictions of time-dependent Kepler (iWeaver et al. 


1978h model predictions (see 


Heger et al. 


20071). The observation-model comparisons indicate 


that the source normally undergoes rapid proto n (rp)-process 


approximately solar composition (i.e. Case 1 of iFujimoto et al 


jumin g of mixed H/He fuel with 


19811. 


Subsequent analysis of a more extensive burst sample showed deviations from the previously 
tight correlation between the flux (measured above 2.5 keV by RXTE) and recurrence time. 
However, simultaneous Chandra and XMM-Newton observations indicated that these deviations 
may result from underestimates of the persistent flux arising from a partial redistribution to lower 
energies, such that the accretion rate-recurrence time relationship remained close to that expected 


^http://burst.sci.monash.edu/wiki/index.php?n=MINBAR.SourceTable 
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theoretically (IThompson et alJl2008h . 


Due to the absenee (so far) of Eddington-limited bursts, the souree distanee has been 


(Barret et al. 

200 

0), while the pea 

c flux of sub-Eddington bursts implies an upper limit of 8 kpe 

(in’t Zand et al. 

1999; 

Kone et al. 

2000 

). By matching the observed burst profiles with Kepler 


numerieal model predietions, 


Heger et al. 


2007h estimated a distane e of (6.07 + 0.18) kpe, 


where €h is the burst emis sion anisotropy faetor. 


Zamfir et al 


data as 


Heger et al. 


(120121) analyzed the same RXTE 


(|2007h to establish mass and radius constraints, as well as an upper limit on 


the distance of 5 .5 x kpe. For these constraints (and also for this paper) the convention of 


Fujimotol (Il988h has been adopted, which defines (and the corresponding value for the persistent 
emission, ^p, which may have a different value) such that the luminosity = An(f^b,pFb,p- 
Thus, ^b,p > 1 implies that emission is preferentially beamed away from the line of sight, so that 
the isotropic luminosity implied from the flux measurements is an underestimate. 

Since its discovery, GS 1826-24 has consistently been observed in a persistent “hard” 
spectral state characterized by a dominant power-law component. Other burst sources are known 
to switch between hard and soft states, the latter associated with higher accretion rates, which 
last for days to months, and are accompanied by changes in burst behavior. Due to the pattern 
described by these so urces in an X-ray color- color diagram, these states are known as the “island” 


and “banana” states (Ivan Paradijs et al. 


19881 see also Galloway et al. 2008). In NS-FMXBs, 


spectral state transitions are thought to involve variations in the accretion flow through a truncated, 
optically thick and geometrically thin disk. In the low-hard (island) state the accretion disk inner 
radius is limited by a hot optically thin quasi-isotropic inner flow, while in the high-soft (banana) 
state, the hot flow v anishes as the disk inner radius extends down to the NS surface, meeting a 


boundary layer (see 


Barret 


2001 


Done et al. 


20071) . These changes in accretion flow geometry 


are related to changes in mass accretion rate, and are thought to affect the burst behavior. As 
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an example, the transient X-ray burster IGR J17473-2721 was observed in outburst in 2008, 
experiencing a remarkable switch from hard to soft state accompanied by a dramatic change 
in burst behavior, w hich demonstrated a h ysteresis in the burst rate as a function of persistent 


bolometric flux (see 


Chenevez et al. 


201 ill) . Another particular effect of the accretion flow on 


the burst behavior is the interaction, during the soft state, of the boundary layer with the NS 
atmosphere, which influences the spe ctral evolution of t he burst emission in a way that is not 


observed during hard state bursts (see 


Kajava et al. 


2014 and references therein). 


On 2014 


(Nakahira et al 


une 8 , GS 1826-24 was detected for the first time in a soft spectral state 


20141 see also Asai et al. 2015), which lasted mo re than two months, a ccording 


to the long-term monitoring by the MAX/ Gas S 


Swift Burst Alert Telescope (BAT; 


Krimm et al. 


it Cam era (GSC; 


Matsuoka et al. 


2009)and the 


20131) . Here we present analysis of NuSTAR 


and Swift target-of-opportunity (ToO) observations of GS 1826-24 triggered in response to this 
unprecedented episode. 


2. Observations and data Analysis 


2.1. MAXI 


The Monitor of All-sky X-ray Image (MAXI; 


Matsuoka et al. 


2009) has been deployed 

I, 


aboard the International Space Station since August 2009. We use publicly available dataj from 
MAX//GSC to examine the long-term 2-20 keV intensity of GS 1826-24. We converted the 


observed GSC count-rate to mCrab units adopting 3.3 + 0.1 count cm ^ s ^ for 1 Cral 


as obtained 


from the average GSC count-rate over the same time interval between October 2013-October 


-http://maxi.riken.jp/top/index.php?cid=l&jname=J1829-237#lc 
"^Equivalent to a flux of (3.2 ± 0.1) x 10“* erg cm“^ s“^ (2-20 keV). 
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2014. 


2.2. Swift 


We utilize daily averaged 15-50 keV intensity measurements 


by BAT 
website ^ 


Barthelmv et al 


Krimm et al. 


2005h on the Swift satellite (IGehrels et al. 


or GS 1826-24, measured 


2004) downloaded from the 


20131) for this analysis. A long-term light eurve was extraeted over the 


same time interval as for the MAXI data. The BAT eount-rate was eonverted to mCrab adopting 1 


Crat 


= 0.22 + 0.008 eountem ^ s 


2 C .-1 


The Swift X-Ray Teleseope (XRT; 


Burrows et al. 


20051), whieh is sensitive to X-ray photons 


in the 0.2-10 keV ban d, observed GS 1826 -24 on 2014 June 20 for 1 ks as a follow-up to the 


report of the soft state (iNakahira et al. 


20141) . On June 24 we requested a longer ToO observation 


with the goal of deteeting X-ray bursts, obtaining an additional exposure of 17 ks. A third 
observation was seheduled to eoineide with out NuSTAR ToO (see below) on June 27, for 1.5 ks. 
All these XRT observations (see Tabled] for details) were exeeuted i n window timing (W T) mode. 


2009h provided 


The raw data were first redueed using the online XRT produets tool (lEvans et al. 
by the Swift team at the University of LeieesteiQ, and with our own analyses, whieh gave eonsistent 
results with the former. Our analyses, whieh are used in the present paper, were performed with 
standard software within HEASOFT v6.16 and CAEDB files from 2014 June 10. 


"http://swift.gsfc.nasa.gov/results/transients/Gingal826-238/ 
^Equivalent to a flux of (1.5 + 0.1) x 10^^ erg em“^ s“^ (15-50 keV). 

'http://www.swift.ac.uk/user_obj ects/ 
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2.3. NuSTAR 


The Nuclear Spectroscopic Telescope Array {NuSTAR', 


Harrison et al. 


2013 ) consists of 


two identical telescopes with a 10 m focal length, focusing X-rays between 3-79 keV using 
depth-graded multi-layer grazing incidence optics. At the focus of the telescopes are Focal Plane 
Modules A and B (FPMA and FPMB), each consisting of a grid of four CZT sensors, with 32 x 
32 pixels. 


NuSTAR performed a ToO observation of GS 1826-24 on 2014 June 27 and 28 for a total 
elapsed time of 108 ks divided into two contiguous datasets with exposures of 13.2 ks and 38.7 ks, 
respectively (Table [U). The NuSTAR data were reduced using the standard NuSTARDAS pipeline 
V 1.4.1 utilizing CALDB files from 2014 October 20. Images obtained from FPMA and FPMB in 
each dataset were used to define source and background extraction regions, both situated on the 
same pixel sensor of the detector. Light curves and spectra of GS 1826-24 were extracted using 
the FTOOLS ‘ ‘nuproducts ’ ’ from a region of 100" radius centered on the source location in 
each module. Based on the NuSTAR point spread function (PSF), this aperture contains 99% of 
the source counts. Another circular region of 120" radius centered about 280" from the source 
was used to measure the sky and instrument backgrounds. The background outside the source 
extraction region is negligible (< 1% of the source counts) below 30 keV. For analysis of the 
persistent emission, we subtracted the full-bandwidth background spectrum, extracted over the 
same time interval as the source spectrum. 
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2.4. INTEGRALUmA-X and RXTEI^CA data from MINBAR 


In this paper we utilize preliminary data from the Multi-INstrument Burst ARehive 


(MIh 


2006 


BAIO- They eonsist o 


and BeppoSAXfWVC (IBoella 


analyses of all bursts deteeted in publie RXTEfPCX (IJahoda et al. 


1997h data through the who 


well as all pubhe dat a from the JEM-X eamera (iLund et al. 


(Winkler et alJ 


2003 


e lifetimes of these missions, as 
I onboard the INTEGRAL satellite 


20031) . through 2014 Deeember. Analysis produets inelude full-range light eurves 


at 0.25 s (1 s) time resolution, fox RX TE (INTEGRAL) , as we ll as time-resolved speetral analyses 


following the proeedures deseribed by 


Galloway et al. 


(2008). 


2.5. Time-resolved spectral analysis 


We extraeted time-resolved speetra eovering eaeh burst observed with Swift and NuSTAR 
and earried out speetroseopy on these data as follows. We first defined time bins using full-energy 
range light eurves at 0.25-s time resolution. We subtraeted the pre-burst level, and defined time 
bins forward and baekward from the time of peak eount-rate sueh that eaeh bin had approximately 
the same number of deteeted eounts. For Swift, the aim was 350 eounts; for NuSTAR, the aim was 
200 eounts eaeh in FPMA/B. The shortest time bin for the Swift burst was 3 s; for the NuSTAR 
bursts, 1 s. Half of the time bins for the NuSTAR bursts were 3 s or shorter. Trial-and-error 
suggests that shorter bins offer no improvement on the speetral fit parameters. 


The NuSTAR burst data w ere signifieantly affee ted by dead time, as is eommonly the ease 


when observing bright objeets (IHarrison et al. 


20131) . This effeet reduees th e deteeted eount-rate 


below that ineident on the deteetors, and so a eorreetion must be applied (see 


Baehetti et al. 


2015t 


The most energetie burst (#3; see seetion reaehed a peak net intensity of approximately 


^http://burst.sci.monash.edu/minbar 
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1200 count which corresponds to almost twice the Crab eount-rate (eorreeted for dead time, 
PSF, and vignetting). At this intensity, and ineluding the pre-burst (persistent) emission, the dead 
time fraetion was about 0.75. At the median eount-rate for all the bins of 200 eount s“', the dead 
time fraction was 0.4. The high dead time fraetion neeessitated the time binning deseribed above 
being performed on the deteeted eounts (rather than the inferred ineident eount-rate). 


We rebinned eaeh speetrum to ensure at least 10 eounts per bin. We fit eaeh speetrum with 
an absorbed blaekbody model, with the neutral absorption fixed at 4 x 10^' em“^ (in’t Zand et al. 
1999). For the Swift speetra, we fit in the range 0.3-10 keV and ineluded a systematie error of 3%, 
as reeommended in the Swift CALDB release note #5|d. For the NuSTAR speetra, we assumed no 
systematic error and fit in the energy range 3-20 keV. 


3. Results 

We show the long term intensity and speetral state history of GS 1826-24 in Fig. [U via 
the joint MAXI and SwiftfQPH light eurves, where the times of the NuSTAR and Swift ToO 
observations are indicated. 

Beginning around MJD 56803 (2014 May 26) the 2-20 keV MAXIjGSC intensity inereased 
over a week-long interval to more than a faetor of two higher than the typieal value of 45 mCrab. 
During this exeursion, the 15-50 keV Swift[QPH intensity was steady. A eloser inspeetion of the 
MAXI light eurve at the orbital resolution reveals that the flare on MJD 56810 was likely due 
to an X-ray burst, and a handful more bursts were detected by MAXI all along the souree soft 
state episode. One of these bursts oeeurred during the time interval eovered by our observations, 
but unfortunately at a time coinciding with one of the NuSTAR orbital data gaps (see ^3.21) : the 
other bursts were separated by more than one day from our Swift and NuSTAR observations. 


^http: //WWW. swift.ac.uk/analysis/xrt/files/SWIFT-XRT-CALDB-09_vl6.pdf 
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GS 1826-24 


o October 2013 - 2014 



Fig. 1.— Daily averaged persistent intensity of GS 1826-24 between 2013 October 1 and 2014 
October 30 as measured by MAXI and 5wz/t/BAT. The data gap between MJD 56640 and 56665 
corresponds to the time the source could not be observed due to instrumental sun-angle constraints. 
The time interval of our SwiftfXKV and NuSTAR observations is indicated by vertical dashed lines 
(MJD 56832 - 56836). Arrows on the time axis indicate the dates of the bursts detected by 
/ATEGiML/JEM-X (see section 14.31) . and the corresponding observation coverage is shown on 
the horizontal 0-line. 
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The MAXI/GSC daily intensity returned to the pre-flare level by MJD 56814, but one day later 
began to increase again, this time accompanied by a steep decrease in the Swift/BAI intensity. On 
MJD 56820 (2014 June 12), the 15-50 keV intensity became below the Swift/BAI detection level, 
and remained below 20 mCrab for the next 30 days. The 2-20 keV intensity was above the typical 
level through to MJD 56850, excluding a 5-day interval beginning MJD 56832 (and coinciding 
with the scheduling of our ToO observations). The 15-50 keV intensity recovered to the typical 


level of approximatelv 
sake of completeness 


10 mCrab over a much more extended period of about 50 days. For the 
we note a previous episode in 2013 October-November during which the 



3.1. Persistent emission 

We investigated the persistent spectrum using 0.3-10 keV SwiftfXKT and 3-78 keV NuSTAR 
spectra. It appears from the long term MAXI light curve shown in Fig. [T]that our Swift and 
NuSTAR observations were performed while the 15-50 keV intensity was still suppressed, but the 
2-20 keV intensity had temporarily returned to « 50 mCrab, roughly consistent with the level 
prior to the flaring activity. 

The source intensity light curve obtained with NuSTAR is shown in Fig. [2l where both the 
variation of the persistent intensity between bursts, and the peak count-rate of six bursts are 
displayed simultaneously. Burst #3 (see below), which has the highest peak intensity, occurred 

^°At the time of writing this paper, we note that a similar soft episode of GS 1826-24 was 
recorded by BAT and MAXI for a duration of about 20 days around 2015 June 3 (MJD 57176), and 
again from 2015 July 9 (MJD 57212) through August 2015. 
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NuSTAR light curve of GS 1826-24 
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Fig. 2.— X-ray intensity of GS 1826-24 measured in the 3-42 keV band by A^M57]4i?/FPMA 
during 2014 June. The persistent emission (at 1000 s) resolution is shown (black symbols and 
histogram, left-hand i/-axis) along with the time of the bursts (dashed red lines); the peak intensity 
is indicated by the length of the lines (right-hand y-axis). The blue arrow indicates the time of the 
MAX/burst (#4). 
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after the longest separation from the previous event (assuming no burst is missed during the 
regular data gaps). The persistent count-rate was steady at approximately 32 count s“' within this 
interval, dropping slightly to a minimum immediately following the burst, and from this point 
rising steadily to a level about 30% higher towards the end of the observation. 


also 


In order to establish a cross-calibration of SwifiPiKI and A^M57]4i?/FPMA and FPMB (see 


Madsen et al 


20151) . we first identified all the times of overlap between the observations 


with the two instruments. There were only two such intervals, between MJD 56835.94105 and 
56835.94730 (duration 540 s), and between MJD 56835.99881 and 56836.00189 (duration 266 s). 
We refer to these two intervals as 01 and 02, respectively. We extracted SwifiPiKT spectra from 
observation 00080751002 over each of these intervals, and NuSTAR/FFMA and FPMB spectra 
from observations 80001005002 (01) and 80001005003 (02). 


We carried out a joint fit of the spectrum for both intervals 01 and 02 simultaneously in 
the range 0. 3-10 keV (XRT) and 3 -40 keV (NuSTAR), with the double Comptonization model 


adopted by 


Thompson et al. 


(|2008l) . We grouped the XRT and NuSTAR spectra to ensure a 


minimum of 10 counts per bin for XRT, and 30 counts per bin for NuSTAR. No source emission 
was detected with NuSTAR above 50 keV. We se t the neutral absorptio n along the line of sight 


with the column density frozen at 4 x 10^^ c m ^ din ’t Zand et al. 


medium abundances (Xspec tbabs model of 


two compTT components in Xspec (lAmaud 


Wilms et al. 


3 


19991) with updated inter-stellar 


2000|) . The composite model consists of 


19961. and references therein), one with a low electron 


temperature kTg and high optical depth r, and the other with a high kTg and low r. The electron 
temperature for the high-kT^ component was effectively unconstrained i n the fits, and so we fro ze 


this value at 20 keV (as measured by RXTE observations in 2002-3; see 


Thompson et al. 


2008). 


The resulting fit, with the spectral parameters tied between the two intervals, gave a reduced 
= 1.011 for 1804 degrees of freedom. The full set of spectral fit parameters are listed in Table 
[2l and the unfolded spectrum and data-to-model ratio for interval 01 are shown in Figure [3l 
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Simpler spectral models, such as a single Comptonization component, do not yield acceptable 
fits for plausible absorption columns. We also included a constant multiplicative factor in order 
to establish any relative flux offset between the SwiftPiKT and NuSTAR instruments. The best-fit 
value of this parameter was 1.022 + 0.015, indicating that the two instruments are consistent 
within their absolute flux calibration. 

As expected based on the SwiftfBAT and MAXI light curves, the spectrum measured by Swift 
and NuSTAR in 2014 June was substantially softer than previous measurements. The electron 
energy for the softer component (^Tg i in Table |2l) was about a factor of two lower, while the 
optical depth ti was similar. Although we cannot constrain the electron temperature for the 
second component, with that parameter fixed at roughly the same value observed previously, the 
optical depth for this, T 2 , was less than half the previous value, indicating a spectrum decreasing 
much more steeply to higher energies. This is illustrated by the comparison with the most 
recent RXTE observation, on MJD 55683.59171, of 10.126 ks duration, with the Proportional 
Counter Units (PCUs) 1,2,4 active (Fig. [3]). At that time the hardness ratio of SwiftfBAT to MAXI 
intensities was > 1.7, compared to the corresponding value in 2014 June of < 0.7. 

We then applied the double Comptonization model to each of the inter-burst intervals for 
bursts #2-8 (see 13.21) . For the interval between bursts 2 and 3, which spans the two NuSTAR 
observations, we extracted for simplicity a spectrum only from observation 8001005003, because 
the average count-rate was about the same; this covers 6.81 hr of the total (13.636 hr) separation. 
We fitted these spectra simultaneously with the double Comptonization model, and experimented 
by trial-and-error, allowing different combinations of parameters to vary between the intervals. We 
first freed each of the Comptonization normalizations, and found that freeing only one additional 
parameter, T 2 , was sufficient to obtain an adequate fit overall, wiihxl = 1.0305 {P = 0.069) for 
4806 DOF (Table O. We used the cflux convolution model component in Xspec to measure the 
unabsorbed model flux within each interval in the 3-25 keV energy range. 
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Fig. 3.— Persistent speetra of Swi/t/XRT and NuSTAR data for GS 1826-24 during interval 01 
of the 2014 June observation, along with the most reeent RXTE speetrum, from 2011 May 2 {top 
panel). The speetra are plotted as (v) to highlight the differenee in speetral hardness. The blaek 
symbols (histogram) show the data (model) for the XRT data, while the red and green symbols 
(histogram) show the data (model) for the NuSTAR FPMA and FPMB data. No souree emission 
was deteeted by NuSTAR above 50 keV. The blue symbols (histogram) show the speetrum observed 
by RXTE on MJD 55683.59171 (obs ID 96306-01-01-03, data only from PCU #2 shown), when 
the souree was still in its hard state. The fitted model for the XRT and NuSTAR data are the best¬ 
fitting double-Comptonization model; for RXTE, a single Comptonization eomponent is used. The 
lower panel shows the data-to-model ratio for the best-fit model with parameters listed in Table |2] 
for the XRT and NuSTAR data. 
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As is customary, we used an “ideal” response to extrapolate the best-fitting speetral model 
outside the in strument bandpass to t he range 0.1-1000 keV, and ad opted this as the bolom etric 


flux (see also 


Thompson et al. 


2008h . Unlike the previous study by 


Thompson et al 


(l2008h . for 


whieh the correction to the bolometric flux based on the absorption was approximately 5%, the 
speetrum during the 2014 June observations was so soft that the eorreetion was cl oser to 35%. We 


estim ate the average unabsorbed bolometrie flux (for comparison to the results of iThompson et al. 
2008) at (2.7 + 0.2) x 10“^ erg em“^ s“^ 


3.2. Thermonuclear bursts 


We deteeted seven type-I (thermonuclear) bursts during our ToO observations in 2014 
June, the first by Swift on MJD 56832.99124, and the remaining six detected by NuSTAR as 
shown in Fig. [2l As mentioned above, one more burst appears in the MAXI orbital light curve at 
MJD 56836.5224 + 405 (M. Serino; private eommunication). We do not have more detailed data 
for this burst, but we note that its time coincides with a NuSTAR orbital data gap. We number 
these bursts #1 to 8 based on their chronological order (see also Table |3]). 


The bursts we observed were significantly shorter than previous bursts, as determined by the 
duration over which the count-rate exeeed 25% of the maximum. The typieal timescales are « 12 s 
(see Table[3]), eompared to 35.9 + 0.4 s for the RXTE bursts in the MINBAR sample (Fig. H]). We 
fitted a one-sided Gaussian to the rising part of each burst and translated the standard deviation to 


the time it takes the Gaussian to rise from 25% to 90% of t 
times the standard deviation; see also iGalloway et al.ll2008 


re peak value (corresponding to 1.206 
. We observed eonsiderable diversity 


among the NuSTAR bursts, both in burst rise time (in the range 2-12 s) and peak intensity (a range 
of a faetor of 8). This ineonsisteney between suceessive bursts is also atypical for this source (e.g. 
Galloway et al. 2004). The brightest NuSTAR burst exhibited the shortest rise time, of 2 s. The 


last burst observed (#8) was also the weakest, and oecurred after the shortest reeurrenee time ever 
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observed in this souree (see below). The rise time for this burst (and also burst #5) was similar to 
the deeay time, so that the burst was almost symmetric in profile. 



Time [s] 


Fig. 4.— Light-curves of the SwifiPiKT burst (#1; 0.2-10 keV) and the six NuSTAR bursts (#2, 3, 
5, 6, 7, 8; 3-42 keV) compared with a burst detected with RXTE (2-60 keV) on 1997 November 
5. The pre-burst average count-rates are subtracted, and the 5w//r/XRT and RXTE burst peaks are 
normalized to the highest NuSTAR peak (burst 3) at about 1200 counts s"'. 

The shortest separation between any observed burst pair was between the last two bursts 
observed by NuSTAR, #7 and 8 in Table[3l These events were observed on MJD 56836.82851 and 
56836.88213, respectively, with a separation of 1.287 hr. Previously, GS 1826-24 has exhibited 
consistently regular bursts, so we tested whether the bursts observed in 2014 June were consistent 
with a regular recurrence time. The separations of the previous three pairs, at 2.1, 3.179 and 
2.082 hr, respectively, are not consistent with the separation for the final pair, nor any integer 
multiple, as expected if bursts were missed in data gaps. However, the final burst observed with 
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NuSTAR was much weaker in both fluence and peak flux (see Table [3]) than the other bursts 
in the same observation, so we consider the possibility that the final burst was the second (or 


third) component of a so-called “short-recu rrence time 


in most sources accreting mixed H/He fuel (IKeek et al. 


mrst” , groups of up to four events seen 


2010t) . Although 1.287 hr is beyond the 


usual range of delays seen for such events, it is possible the last burst followed more closely 
another event which fell in the data gap which ended just 7.81 min earlier. In that case, we should 
discount the final burst, and consider only the four previous ones (including the MAXI burst). 
The separation between the successive pairs of bursts were related in approximately a 3:2 ratio, 
suggesting that the bursts were occurring regularly every 1.05 hr. However, if that were the case, 
the expected time of one of the missing bursts between the observed events #5 and 6 fell in the 
middle of an observation interval in which no bursts were observed. Thus, we can rule out regular 
bursting during the time interval covered by our observations at high confidence. 


3.3. Burst energetics and spectral variations 


We carried out time-resolved spectral analysis as described in ^2.51 We found an adequate fit 
to each time-resolved net burst spectrum (with the pre-burst emission subtracted as background) 
using an absorbed blackbody model. The resulting distribution of reduced values is shown in 
Fig. [51 The maximum value for any of the fits was 1.34; this is consistent with expectations for a 
good fit given the number of degrees of freedom. We list the burst spectral parameters in Table [3l 


Figure 6 shows time-resolved spectroscopic results for the Swift and NuSTAR bursts. The 
time-resolved spectroscopic analysis of the brightest burst, #3, indicates the characteristic 
evolution of a photospheric radius expansion (PRE) burst, with a local maximum observed in the 
blackbody normalization at the same time as a minimum in the blackbody temperature (Fig. Ibbl) . 
he presence of P RE in other sources is strongly correlated with the source being in a soft state 


(Muno et al. 


20041) . 
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Fig. 5.— Distribution of reduced-;^f^ fit statistic for 76 spectra covering the six bursts observed by 
NuSTAR from GS 1826-24. The smooth curve overplotted is the expected distribution assuming 
the adopted model (an absorbed blackbody with fixed neutral column density) is correct. The 
smooth curve is calculated for the average number of degrees of freedom in the fit (50). 
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The peak flux reached during this burst was (40 + 3) x 10“^ erg cm“^ s“^ This value is a 
factor of 1.42 higher than the mean peak flux of the non-PRE bursts observed by RXTE since 
2000, of (28.4 + 1.2) x 10“^ erg cm“^ s“\ and a factor of 1.24 higher than the peak flux of the 
brightest burst yet observed from the source: (32.6 + 1.0) x 10“^ erg cm“^ s“^as qutoed from the 
MINBAR database. 



0 10 20 30 40 


Time from burst start (s) 

Fig. 6a.— Time-resolved spectroscopy of burst #1, observed by SwiftPiRT on MJD 56832.99124. 
The top panel shows the inferred bolometric luminosity, assuming a distance of 5.7 kpc. The mid¬ 
dle panel shows the best-fit blackbody temperature, and the lower panel shows the normalization. 

We integrated over the measured fluxes to give the fluence Eh for each burst, and computed 
the burst timescale r as the ratio of the fluence to peak flux, i.e. r = E^/Fpeak- A measure of 
burst energetics is given by a = A4ecFpers/Ffo, the ratio of persistent and burst fluences (e.g. 
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0 10 20 30 40 10 20 30 40 10 20 30 40 


Time from burst start (s) 

Fig. 6b.— Same as Fig. 6a for bursts #2, #3, and #5, observed by NuSTAR. Note the moderately 
strong radius expansion of burst #3 during the first four seeonds of the rise, eoupled with a deerease 
in kTbb- 


Lewin & Joss 


1983h . In this expression for a, Fpsrs and At„ 


are the average bolometrie persistent 


flux and the waiting time sinee the last preeeding burst, respeetively. We use the MAXI burst so as 
to better eonstrain the or-value of the following NuSTAR burst (#5), although this is obtained with 
a relatively high uneertainty due to the approximate knowledge of the time of the MAXI burst. 
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0 10 20 30 40 10 20 30 40 10 20 30 40 


Time from burst start (s) 

Fig. 6c.— Same as Fig. 6b for bursts #6, #7, and #8, observed by NuSTAR. 

4. Discussion 

The 2014 June soft spectral state of GS 1826-24 was the first ever recorded for this 
well-studied source, and it revealed a number of new observational features, including the first 
burst exhibiting photospheric radius expansion, and weak, irregular bursting behavior, including 
the shortest burst interval (1.29 hr) measured to date. 
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4.1. The source distance 


The brightest burst observed with NuSTAR, #3, exhibited speetral evolution eonsistent with 
PRE, thought to indieate the burst flux reaehing the Eddington limit. Assuming the peak flux 
eorresponds to the Eddington luminosity for an atmosphere with solar eomposition, and taking 
into aee ount the effeets of grav itational redshift at the surfaee of a 1.4 M©, 10 km radius neutron 


star (e.g. 


Galloway et al. 


2008h the inferred distanee is (5.7 + 0.2) A kpe, where represents the 


possible anisotropy of the burst emission (see ^4.31) . In at least one other syste m, 4U 1636-53 6, 


whieh is thoug 


It to a eerete mixed H/He fuel as assumed for GS 1826-24 (see 


Galloway et al. 


Bildsten 


20001 : 


20041) . the effeetive E ddington limit is 


appropriate for a pure He atmosphere (iGallowav et al. 


implied distanee is (7.4 + 0.3) kpe. Eurther, 


thoug ht to be instead the higher limit 


20061). If we ins tead adopt that value, the 


Kuulkers et al. 


(120031) measured the Eddington 


luminosity for a group of EMXBs with independently known distanees from their globular eluster 
host as (3.79 + 0.15) x 10^^ erg s"'. Based on this value the implied distanee is (8.9 + 0.4) 
kpe. 

These larger distanees are problematie for several reasons. Eirst, the non-PRE bursts observed 
previously reaeh an average maximum flux only a faetor of 1.42 lower than burst #3, implying 
that the non-PRE bursts exeeed the Eddington limit for mixed H/He fuel. This also seems to 
be th e ease for 4U 1636-5 36, whieh infrequently shows PRE bursts eonsistent with the H/He 


limit (Galloway et al. 


2006h . However, for GS 1826-24 the He-derived distanees also exeeed the 


upper limit of 5.5 x A kpe derived frorn eomp aring the non-PRE burst light eurves to Kepler 


numerieal model predietions (IZamfir et al. 


2 OI 2 I) . Thus, though we eannot absolutely rule out 


other possibilities, we adopt a distanee o f (5.7 + 0.2) kpe as this is the only one that satisfies 


the eonstraints obtained by 


Zamfir et al. 


(120121) . and we eonelude that the effeetive Eddington 


limit for GS 1826-24 is for mixed H/He fuel. 
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4.2. The persistent spectral state 


Since previous observations of GS 1826-24 have eonsistently found the souree in the hard 
(island) speetral state, we diseuss here to what extent the 2014 June observation is distinet from 
that state. Extensive previous RXTE observations of other “atoll” elass LMXBs (so named 
beeause of their eharaeteristie pattern in X-ray eolor-eolor diagrams) find that the hard and soft 
X-ray eolors (defined as the ratio of eounts between pairs of energy bands - for RXTE, the energy 
bands used were 8.6-18.0 and 5.0-8.6 keV for the hard eolor, and 3.6 -5.0 and 2.2-3.6 keV for 
the soft eolor) of these sourees define an are or a Z-shaped traek (e.g. 


Galloway et al. 


2008). 


Unfortunately, beeause GS 1826-24 was never observed by RXTE to go into a soft state, it’s 
eolor-eolor diagram is not well defined, and instead all observations eluster around a soft eolor 
value of 1.69 + 0.04, and hard eolor 0 . 865 + 0.013. These values, extraeted from the eatalog of 


RXTE observations of 


Galloway et al. 


(l2008h . are eorreeted for the PC A gain, whieh varied over 


the mission; the eorresponding values prior to the gain eorreetion for the epoeh elosest to the end 
of the mission would be 15% lower in soft eolor, and 4% higher in hard. 

We estimated the eorresponding eolors for the 2014 June observation in the same energy 
bands used for the RXTE analysis, and for the most reeent gain epoeh. We ereated a simulated 
persistent speetrum in Xspec, adopting the best-fit speetral model with a response ealeulated 
for a late-epoeh RXTE observation. The estimated PCA eolors for the souree in 2014 June are 
1.189 (0.476) for soft (hard) eolor. In other atoll sourees, a signifieant deerease in both soft and 
hard eolor is assoeiated with a transition to the “banana” or soft speetral state. Although it is 
impossible to be eertain in the absenee of a well populated eolor-eolor diagram for GS 1826-24, 
the speetral measurements strongly support a state transition similar to that seen in other atoll 
sourees. Furthermore, although the higher aeeretion rate that might be implied by the speetral 
state transition is not supported by the estimate of t he bolometrie flux, sue h diserepaneies are also 


well-known in other atoll sourees (see e.g. Fig. 6 of 


Galloway et al. 


2008). 
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For other atoll sources, the soft “banana” persistent spectral state is usually interpreted as 
indicating a higher accretion rate, and naively, the higher average burst rate for GS 1826-24 
during 2014 June would seem to support this interpretation. However, the inferred bolometric 
persistent flux level of (2.7 + 0.2) x 10^^ erg s“^ is in fact in the middle of the range of 


bolometric fluxes that the source has been observed at historically (e.g. 


Thompson et al. 


2008h. 


Thus, we find no evidence to support a markedly different accretion rate, unless the radiative 
efficiency (or perhaps the persistent emission anisotropy) has changed markedly. The Swift and 
NuSTAR observations fell between two much higher peaks of the MAXI light curve (Fig. [B; it 
seems likely that the source could have been up to a factor of two brighter still in the soft state, 
just a few days before or after. 


We note that 


Ji et al. 


(120141) report the diminution of the hard X-ray persistent emission 


during GS 1826-24 bursts observed by RXTE. These authors explain such hare 
as due to the cooling of the hot corona by the soft X-ray burst photons (see also 


X-rav shortages 

Ji et al. 

2015 

). We 


tested for similar variation in the hard X-ray emission during the six bursts observed by NuSTAR, 


but the source intensity above 30 keV was persistently so weak (on 
could not find any significant variation. This may be consistent with 


V a few count s/s) that we 


Ji et al. 


(120141) results, as our 


observations occurred when GS 1826-24 was in a soft state during which a negligible corona or 
hot accretion flow is supposed to be present. 


4.3. The bursting regime 


GS 1826-24 has so far been characterized by consistently regular > 100 s long bursts 


recurring at approximately periodic intervals betwe en 3.56 and 5.7 


almost linear function of the source persistent flux (IGalloway et al. 


hr, v arying inversely as an 


20041). Apart from the burst 


detected by MAXI on MJD 56810 (see ^3]), we found no observations in the few weeks preceding 
or succeeding the soft episode. However, JEM-X detected two bursts before the soft spectral 
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episode, in 2013 October 23 and 28 (MJD 56588 and 56593), and afterwards, on 2014 September 
15 and Oetober 25 (MJD 56915 and 56955), respeetively. As shown in Fig. [H the four JEM-X 
bursts oeeurred while the souree speetral state was rather stable, with the (BAT / MAXJ) hardness 
ratio eonsistently about 1.6 (+10%). All four bursts before and after 2014 June show similar 
shapes and durations as other bursts from GS 1826-24 previously observed with JEM-X in the 
3-25 keV energy range, i.e. longer rise times (from 6 to 9 s) and timeseales (between 40 and 
60 s). While it is not possible to infer the burst rate with sueh widely separated burst deteetions, 
the long burst timeseales and eharaeteristieally hard persistent speetral state strongly suggests that 
GS 1826-24 was exhibiting its normal burst behavior up to 2014 June, and following. 


Based on previous measurements (IThompson et al. 


20081) . at the flux level seen in 2014 June 


we would expeet regular, eonsistent bursts at a recurrenee time of At « 4 hr, and with a ~ 30-40. 
Instead, we found mueh weaker, inhomogeneous bursts, with fluenees at most one third of the 
typieal value measured in the past, and eorrespondingly higher or-values. Given the laek of 
regularity in the bursting, and the presenee of gaps between eaeh of the burst pairs, the measured 
cr-values must be eonsidered upper limits only, so we eannot rule out lower values, eonsistent with 
the usual mixed H/He burning. However, we ean determine a lower limit on the amount of H in 
the burst, based on the assumption that all the aeereted fuel is burned during the burst: 


a = 58 


M 


Qn 


1.4A/q/\ 10km/ \4.4MeVnueleon ' 




1.31 


-1 


( 1 ) 


(note that the expression in 


Galloway et al. 


2008 


omits the redshift faetor) where M, R are the 


mass and radius of the neutron star, 1 + z = (1 - is the surfaee gravitational redshift, 

and 2nuc = 1.6 + 4X MeV nueleon"', where X is the hydrogen fraetion averaged over the burning 
layer. 


The ratio between anisotropy for the pe rsistent and bur st emissions that appears in Eg. [J 


has been estimated as 1.55 for GS 1826-24 (iHeger et al. 


2007h . The modeling of lEujimotol (119881) 


suggests that a system inelination of 75° is required to give this value of the relative anisotropies. 























-29- 


and further implies that = 1.32 and = 2.04. These values indieate that both the burst and 
persistent flux are preferentially beamed away from our line of sight, and the inferred isotropie 
luminosities will underestimate the true values. The eombined effeet on the measured cr-values 
will be to underestimate the true value by a faetor of 1.55. 


Due to the data gaps falling between eaeh pair of bursts detected by NuSTAR, each of the a 
measurements in Table[3]is an upper limit on the true value. The most constraining value should be 
the minimum, which is obtained for burst #5, although with relatively large uncertainty due to the 
lack of absolutely exact timing of the MAX/burst. The estimated cr = 54 + 8, with the uncertainty 


dominated by the separation from the MAX/burst, of 2.1 + 0.3 hr. The corresponc 
Qnuc, a ssuming the range of gravitational redshift 1 + z = 1.19-1.28 estimated by 


ing value of 


Zamfir et al. 


(120121) . would be in the range 2.8-3.0 MeV nucleon \ implying in turn a hydrogen mass fraction 


at ignition of X = 0.3-0.35. We further note that, while low values of X may arise from steady 
hydrogen burning prior to the burst, there has been insufficient time to reduce it to this degree. 
Indeed, assuming solar CNO metallicity Zcno^ the time to burn all the hydrogen at the base of the 
layer is 

4x = 1 1{Zcno/0.02)-\Xo/OJ) hr (2) 


where Xq is the accreted H-fraction (e.g 


Galloway et al. 


20041) . For solar accreted composition. 


there is insufficient time between bursts to reduce the average H-fraction in the fuel layer to 
explain the or-value of burst #5. 


Another way to understand the discrepancy is by considering the column depth of material 
ignited during each burst, given by: 


Ub 


AncfEb 
3.0 X 10^ 


(1+z) 


lO'^ergcm^^ 

Qnuc 


4.4 MeV nucleon" 



(3) 


( 4 ) 
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again assuming implicitly that all the accreted fuel is burned in the burst. We set the value of Q^uc 
based on the assumed average fuel composition X resulting from hot-CNO burning between the 
bursts, i.e. X = Xo(l - 0.5At/tex) (with the factor 0.5 arising because the burning takes place at the 
base, and we average X over the entire column). We compare this with the column depth accreted 
between two bursts separated by a time interval At, which is AM = mAt, where m = M/AnR^ 
is the mass accretion rate per surface area on the NS. Given the estimated bolometric persistent 
flux Fpers of (2.4-2.9) X 10“^ erg cm^^ s“', and the inferred anisotropy factor, we estimate the 
accretion rate at 12-15% of the Eddington rat e (at a distance of 5.7 kpc). This accretion rate 


is in the range expected for Case 1 burning of 


Fujimoto et al.l (Il98lll . corresponding to mixed 


H/He burning triggered by thermally unstable helium ignition. Assuming that the nuclear burning 
is completely conservative, one would expect yb to be close to AM. Instead, we find that yb 
consistently underestimates AM, even for the relatively close pairs of bursts (Fig. 7). 


We infer that the assumption of conservative burning (i.e. that all accreted fuel is burned 
during the bursts) cannot be reconciled with the data, implying that some other process is reducing 
the available fuel prior to ignition. This burning appears to preferentially removing hydrogen, 
based on the short burst timescales and since the reduction is in excess of the normal steady 
hot-CNO burning. 


We conclude that GS 1826-24 bursts inefficiently in the soft state, and igniting fuel with 
significantly lower hydrogen fractions than the previously inferred solar value. The burst intervals 
are too short for the lower H-fractions to arise purely by hot-CNO burning between the bursts, 
unless the CNO metallicity is of order ten times solar. Steady burning of accreted fuel in addition 
to hot-CNO burning would explain both the relative inefficiency of the thermonuclear bursts, and 
would also provide an extra fuel source to explain the relatively low ignition columns. 


The results of our spectral analysis of the persistent emission compared to previous 
observations indicate a softening, but at a similar inferred accretion rate. Such a softening of the 
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Fig. 7.— Comparison of the estimated accreted column AM, calculated from the inferred accretion 
rate and the inter-burst interval, with the ignition column yt, based on the measured burst fluence 
and the assumed Qnuc arising from the effects of hot-CNO burning between the bursts. Both param¬ 
eters are measured in the local (neutron star) frame, assuming 1 -l- z = 1.28. The dashed line is the 
line of equality, at which point the burning would be conservative. Note that the bursts consistently 
ignite at columns well below that accreted. 
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spectrum would normally be explained by a transition from the usual truncated accretion disk, 
with optically thin inner flow, to an optically th ick flow passing throu gh a boundary layer, as is 


commonly observed in other LMXBs (see, e.g.. 


Barret & Olive 


20021) ■ However, for GS 1826-24 


in 2014 June, this transition is not supported by the data, since the optical depths ti, T 2 for both 
components are lower than in the hard state. Some caution is required in interpreting these 
parameters alone, as they are strongly anticorrelated with the corresponding electron temperatures 
kTe i, kTe 2 , and we fix the latter at 20 keV. Since the evidence for increased mass accretion rate 
is weak, we further attribute the markedly different burst behavior, also to the change in disk 
geometry. Although it is presently not understood precisely how the disk geometry can affect 
the burst behavior, the manif estly different burst pro perties in the soft state has been observed 


in several other sources (e.g. 


Comelisse et al. 


explored in the literatur e (see, e.g.. 


the hard state (see, e.g.. 


Womel et al. 


in’t Zand et al. 


2011 


20031). and this interaction i s increasingly being 


2013 


Kajava et al. 


2014 


Ji et al. 


2015|), also in 


20131). 
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Table 1. Log of Swift and NuSTAR observations of GS 1826-24 in 2014 June 


Date 

MJD 

Instr. 

Obs. ID 

Time range 

(UT) 

Exposure 

(ks) 

No. 

bursts 

2014 June 20 

56828 

SwiftfKm 

00035342005 

18:53-19:10 

0.982 


2014 June 24 

56832 

SwiftfKm 

00035342006 

14:06-00:19" 

16.65 

1'’ 

2014 June 27 

56835 

NuSTAR 

80001005002 

15:36-23:30 

13.2 

1 



SwiftfKm 

00080751002 

22:34-00:13" 

1.501 




NuSTAR 

80001005003 

23:30-22:30" 

38.7 

5 


^End time is on the following day. 


'’The peak count-rate of the SwiftfKKV burst is - 155 count s 







Table 2. Persistent spectral fit parameters for GS 1826-24 in 2014 June 


Double Comptonization model 
Parameter 


Units 


01 &02 


t2-t3 


ti-tA 


Interval^ 

tA-h 






ti-h 


Nh 

kTa,i 

kT,,i 

ri 

kTo,2 

kT,,2 

T2 

xl (DOF) 

Absorbed flux (3-25 keV) 
Unabsorbed flux’’ (0.1-1000 keV) 


lO^’ cm-2 
keV 
keV 

keV 

keV 


10 ® erg cm ^ s ' 
10“® erg cm“^ s“' 


0 092+*’’”^ 

u.u?z Q015 

3 04+’’ ’* 

- 0.16 

4 7Q+0.I8 
+■'+- 0.17 

0 404+*’’” 6 
+’■+”+- 0.015 


0.79 ± 0.05 0.390+“|6 

1.011 (1804) - 

1.1646 ±0.0017 
2.539 ± 0.004 


- (4.0)- 

_0 770+0.018_ 

”■■’'+- 0.024 

-4.2 ± 0.2- 

-2.65 ±0.14- 

- 1.592 ± 0.01- 

- (20)- 

0.494 ±0.016 0.372 ±0.018 0.083 ±0.017 0.064+“[* 0.037+“[* 

- 1.0181 (5650)- 

1.0727 ±0.0016 1.100 ±0.003 1.178 ±0.002 1.183 ±0.003 1.236 ±0.005 
2.283 ±0.003 2.381 ±0.006 2.653 ± 0.005 2.689 ± 0.006 2.782 ±0.011 


’’Time interval between bursts #i and #1+1. 

’’Extrapolated, assuming an ideal response. The flux is calculated as the mean of the fluxes for the models over each of the Swift/XKT and NuSTAR spectra, and the uncertainty is calculated 


as the standard deviation. 
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Table 3. Properties of thermonuelear bursts from GS 1826-24 deteeted in Swift and NuSTAR 

observations in 2014 June 


Burst 



Start time 

At 

Rise time^* 

Timescale^ 




no. 

Instr. 

Obs. ID 

(MID) 

(hr) 

(S) 

(s) 

Peak flux’’ 

Fluence‘^ 

Q,d 

1 

Swiyr/XRT 

00035342006 

56832.99124 


1.9 ±0.3 

11.3± 1.6 

40+80 

^^-20 

0.21 ±0.08 


2 

NuSTAR 

80001005002 

56835.69484 

64.89 

2.4 ±0.2 

12.3 ± 1.1 

13.8 ± 1.1 

0.187 ±0.006 


3 

NuSTAR 

80001005003 

56836.26299 

13.636 

1.13 ±0.08 

8.6 ±0.6 

40 ±3 

0.370 ± 0.009 

337 ±8 

4 

MAXI 


56836.5224'= 

6.2 ± 0.3 






5 

NuSTAR 

80001005003 

56836.60928 

2.1 ±0.3 

5.52 ±0.19 

11.8 ±0.9 

25.1 ± 1.9 

0.335 ± 0.008 

54 ±8 

6 

NuSTAR 

80001005003 

56836.74176 

3.179 

3.35 ±0.16 

12.2 ±0.9 

27 ±2 

0.352 ± 0.009 

86 ±2 

7 

NuSTAR 

80001005003 

56836.82851 

2.082 

3.42 ±0.17 

12.6 ± 1.1 

21.8 ± 1.6 

0.280 ± 0.007 

72.1 ± 1.8 

8 

NuSTAR 

80001005003 

56836.88213 

1.287 

3.4 ±0.3 

12.4 ± 1.4 

6.6 ± 0.6 

0.083 ± 0.005 

156 ±9 


^Measured from the count-rate burst light curve from 25% to 90% of the peak value in the relevant energy band. 
Extrapolated peak bolometric flux in units of 10“^ erg cm“^ s“^. 

‘^Integrated burst bolometric fluence in units of 10“^ erg cm“^. 

‘^As every burst interval was interrupted by at least one data gap, the a-values must formally be considered upper limits. 
®This burst is recorded from the public MAXI orbital light curve. 
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